Fractional charges and quantum phase transitions in imbalanced bilayer quantum Hall 

systems 
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We extend the Composite Boson theory to study shghtly im-balanced bi-layer Quantum Hall sys- 
tems. In the global 17(1) symmetry breaking excitonic superfluid side, as the imbalance increases, the 
system supports continuously changing fractional charges. In the translational symmetry breaking 
pseudo-spin density wave (PSDW) side, there are two quantum phase transitions from the com- 
mensurate PSDW to an in-commensurate PSDW and then to the excitonic superfluid state. We 
compare our theory with experimental data and also the previous microscopic calculations. 



Spin-polarized Bilayer Quantum Hall systems at total 
filling factor vt = ^ have been under enormous experi- 
mental and theoretical investigations over the last decade 
[1]. When the interlayer separation d is sufficiently large, 
the bilayer system decouples into two separate compress- 
ible V = 1/2 layers [2]. However, when d is smaller than a 
critical distance dci , even in the absence of interlayer tun- 
neling, the system undergoes a quantum phase transition 
into a novel spontaneous interlayer coherent incompress- 
ible phase which is an excitonic superfluid state in the 
pseudospin channel. [3-5]. 

Although the ESF phase and FL phase at the two ex- 
treme distances are well established, the picture of how 
the ESF phase evolves into the two weakly-coupled FL 
states as distance changes is still not well established. 
The experiment [6] discovered that although there are 
very little dissipations in both the ESF and FL, there 
is a strong enhancement of drag and dissipations in an 
intermediate distance regime. These experimental obser- 
vations suggest that there maybe an intermediate phase 
separating the two phases. If there is indeed such an 
intermediate phase, even it exists, the nature of the in- 
termediate phase remain very important open problems. 
Using Hartree-Fock or trial wavefunctions approxima- 
tion, many authors [7] proposed different kinds of trans- 
lational symmetry breaking ground states as candidates 
of the intermediate state . In a recent preprint [9], I 
constructed an effective action to study the instability 
driven by magneto-roton minimum collapsing at a finite 
wavevector in the pseudo-spin channel. I showed that the 
instability leads to a pseudo-spin density wave (PSDW) 
with a square lattice structure at some intermediate dis- 
tances. The theory puts the ESF state and the PSDW 
state on the same footing and describe the universality 
class of the quantum phase transition between the two 
states. In the presence of disorders, the properties of the 
PSDW are consistent with all the experimental observa- 
tions [6] in the intermediate distances. 

In this paper, we use the CB theory developed in [9] to 
study the effects of imbalance in BLQH. This theory puts 
spin and charge degree freedoms in the same footing and 
can also be used to explicitly bring out the spin-charge 



connection and classify all the possible excitations in a 
systematic way. As shown in [9], there are two critical 
distances and three phases in balanced case: < d < dci, 
the system is in the ESF, dci < d < dc2, it is in the PSDW 
with square lattice structure, d > dc2, it melts into two 
weakly coupled FL. We find that in the ESF side, the im- 
balance is irrelevant, but the merons carry continuously 
changing fractional charges. While in the PSDW side, we 
map the square lattice PSDW at the balanced case into a 
hard core bosons hopping on a square lattice at half filling 
with long range interactions, then adding a small imbal- 
ance in the BLQH corresponds to adding a small chemical 
potential in the boson model. Through this mapping, we 
find the imbalance drives two quantum phase transitions: 
the first one is from the commensurate PSDW state to 
an incommensurate PSDW (IC-PSDW), the second one 
is from IC-PSDW state to the ESF state. Both transi- 
tions are first order transitions. We compare our phase 
diagram with the recent experimental data [13]. We also 
discuss briefly the effects of disorders and compare our 
results with previous results achieved from microscopic 
calculations. 

Consider a bi-layer system with Ni { N2 ) electrons in 
top ( bottom ) layer and with interlayer distance d in the 
presence of magnetic field B — V x A: 



H — Hq + Hint 

Hq = d'^xcl^{x 



■Ca {x) 



2m 



Hint ^ - J d xd x'5pa{x)Vai3{x - x')Spi3{x') (1) 

where electrons have bare mass m and carry charge — e, 
Ca,a — 1,2 are electron operators in top and bottom 
layers, Spa{x) — cl^{x)ca{x) — Pa,ct — 1,2 are normal 
ordered electron densities on each layer. The intralayer 
interactions are Vn — V22 — e^/er, while interlayer in- 
teraction is V12 = V21 = e'^ /eyr^ + d? where e is the 
dielectric constant. 

Performing a singular gauge transformation 4>a{x) — 

J d X arg(x-x )p{x ) ^^^^-^ whcrC p{x) = Ci{x)ci{x) + 

C2{x)c2{x) is the total density of the bi-layer system. We 
can transform the Hamiltonian Eqn.l into a Lagrangian 
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of the Composite Boson (f>a coupled to a Chern-Simon 
gauge field [8]. We can write the two bosons in terms 
of magnitude and phase (j)a = \/pa + ^Pa^^^" ■, then after 
absorbing the external gauge potential A into a, we get 
the Lagrangian in the Coulomb gauge: 

C = i5p+{\dr6+ - ao) + ^[^V^+ + - V2)V6. 



where m± = mi ± m2. 

We can classify all the possible topological excitations 
in terms of (g, m_) in the following table: 



{mi, 1712) 


(1,0) 


(-1,0) 


(0,1) 


(0,-1) 


rri- 


1 


-1 


-1 


1 




1 


-1 


1 


-1 


q 






V2 


-1/2 



27r 



ao(V X a) 



+ '-6p.dr0- + ^(^V^_)2 + ^5p.V.{q)5p. - hjp. 

where 5p+ = dpi + 5p2,Sp^ — — p2,9± = 9i ± 62, 
they satisfy commutation relations [Spa{x),9(-){x')] = 
2ih6ai3S{x — x'),a,P = ±. p = pi + P2 and / — 'iviV2 
which is equal to 1 at the balanced case, V± = Yii^Yis.^ 
hz = V-p- = V- {pi — P2) plays a role of a Zeeman field. 

(1) OfF-dmgonal algebraic order and Spin-wave exci- 
tation: Neglecting vortex excitations, dropping a linear 
derivative term in 9- which is irrelevant in the ESF state 
and expanding the second term in Eqn.2 which includes 
the coupling between the spin sector and the charge sec- 
tor, we can get the complete forms of the three propaga- 
tors [8]: < 9+9+ >, < 9-9- >, < 9+9- >. 

Performing the frequency integral of the + propaga- 
tor carefully, we find the leading term of the equal time 
correlator of 9+ stays the same as the balanced case 
< 9+{-q)9+{q) >= 2 x |j-|-0(i) which leads to the same 
off-diagonal algebraic order exponent 2 as in the balanced 
case. In the 9,0; —> limit, we can extract the lead- 
ing terms in the 6-9- propagator < 9-9- >= ^2^^^2^^2 
where v is the spin wave velocity in the balanced case and 
we can identify the spin wave velocity in the im-balanced 
case vf^ = fv^ = Avii'2V^ = ^yi{l — vijv^ which shows 
that the spin-wave velocity attains its maximum at the 
balanced case and decreases parabolically as im-balance 
increases. 

As analyzed in detail in [8], it is hard to incorporate 
the Lowest Landau Level (LLL) projection in the CB ap- 
proach, so the spin wave velocity can only be taken as a 
phenomenological parameter to be fitted to experiments. 
But as shown explicitly in the next subsection, the CB 
theory can be used to classify topological excitations cor- 
rectly and systematically . 

(2 ) Topological excitations: There are following 4 kinds 
of topological excitations: A9i = ±27r, A92 = or 
A^i = O,A02 = ±27r. Namely (mi,TO2) = (±1,0) or 
(mi,m2) = (0, ±1) where mi and m2 are the winding 
numbers in 9i and 6*2 respectively. The fractional charges 
can be determined from the constraint V x a = 27r5p and 
the finiteness of the energy in the charge sector: 

g = J- / a • dl = X i (f[V0+ -h {ui - 1/2)^9-] ■ dl 
2-K J 2-K 2 I 



Tabic 1: The fractional charge in im- balanced case 

(2) The finiteness of the energy in the spin sector dictates 
that any finite energy excitation must have m_ = 0, so 
the merons listed in table 1 are confined into the follow- 
ing two possible pairs at low temperature. (1) Charge 
neutral pairs: (±z/i,±l) or (±1/2, Tl)- The NGM will 
turn into these charge neutral pairs at large wavevec- 
tors. The pairs behave as bosons. (2) Charge 1 pair 
{vi, 1) -I- (z^2,-l) or charge -1 pair {-Vi, -1) -|- {-V2, !)• 
The pairs behave as fcrmions. They maybe the lowest 
charged excitations in BLQH and the main dissipation 
sources for the charge transports. 

(3) Im-balancc driven quantum phase transitions in 
the PSDW side: As said in (1), the imbalance is irrele- 
vant in the ESF side, but we expect it is important in the 
PSDW side. Starting from the ESF side, it would be use- 
ful to calculate how the mageto-roton minimum depends 
on the imbalance [12,8], unfortunately this is beyond the 
scope of the CB theory. Here we take a different strategy: 
starting from the PSDW side and studying how a small 
imbalance affects the PSDW. If the imbalance is suffi- 
ciently small, we expect the C-PSDW at the balanced 
case in Fig. 2a is a very good reference state. Because it 
is a square lattice with the "up" pseudo-spins taking sub- 
lattice A and the "down" pseudo-spins taking sublatticc 
B, it is reasonable to start from a lattice model from the 
PSDW side. If we think the PSDW as a charge density 
wave ( CDW ) of bosons at half filling on a square lattice, 
then we can view the ESF to the PSDW as a superfiuid 
to CDW transition in a boson Hubbard model of hard 
core bosons near half filling hopping on square lattice of 
bosons: 

H = -tY^ {blbj + h.c.) - p'Y^Ui 

<ij> i 

+ Vi ^ (n, - l/2)(nj - 1/2) + 1/2 ^ 71^71^ + ••• (4) 

«ik» 



where S'^ = = c\c2,S = b = clci,S^ = ^{c\ci — 
C2C2) = b^b — 1/2 are the pseudo-spin density and bo- 
son operators. At the total filling factor = 1, we can 
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1. 77,: 



= -[771+ -I- (l/i - i/2)m_] 



(3) 



impose the local constraint c[ci 
is the boson density, t is the nearest neighbor hopping 
amplitude, Vi , V2 are the nearest and next nearest neigh- 
bor repulsive interactions between the bosons. The • • • 
may include further neighbor interactions. Because of 
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the long-rang Coulomb interaction in Eqn.2, it is impor- 
tant to keep all the long-range interactions in the lattice 
model Eqn.4. If the chemical potential 1-1 = 0, the bosons 
are at the half filling < rii >= 1/2 which corresponds to 
the balanced case v = 1/2. The particle- hole symmetry 
of Eqn.4 corresponds to the Z2 exchange symmetry of 
the BLQH. If the chemical potential /./, ^ 0, the bosons 
are slightly away from the half filling which corresponds 
to the slightly imbalanced case. 

The boson Hubbard model Eqn. 4 in square lattice 
at generic commensurate filling factors f = p/q { p,q 
are relative prime numbers ) were systematically studied 
in [10] by performing the charge-vortex duality trans- 
formation. Recently, we applied the dual approach to 
study reentrant supersolids and quantum phase transi- 
tions from solids to the reentrant supersolids on extended 
boson Hubbard models Eqn.4 at in-commensurate filling 
factors in bipartite lattices such as square and honeycomb 
lattices. In the following, we apply the results achieved 
in [11] to the present problem on square lattice at and 
slightly away from half filling. At q = 2, there are two 
dual vortex fields ijja and -06 • Moving slightly away from 
half filling / = 1/2 corresponds to adding a small mean 
dual magnetic field i?~5/ = / — l/2in the dual ac- 
tion. The most general action invariant under all the 
MSG transformation laws upto quartic terms is [10,11]: 



where pA = ipa^Ja should be interpreted as the vacancy 
number, while the vortices in its phase winding are inter- 
preted as boson number. Of course, a negative imbalance 
can simply achieved by a particle hole transformation 
V't -5f in Eqn.6. 

Eqn. 6 has the structure identical to the conventional 
q = I component Ginzburg-Landau model for a type II 
" superconductor " in a "magnetic" field. It was well 
known that as the magnetic field increases, there are two 
first order phase transitions: H < Hd, the system is in 
the Messiner phase, Hd < H < Hc2, it is in the vortex 
lattice phase, H > Hc2 it is in the normal phase. In the 
present boson problem with the nearest neighbor inter- 
action Vi > and further neighbor interactions in Eqn.4 
which stabilizes the CDW state at / = 1/2 ( Fig. 2a ) and 
IC-CDW state at / = l/2+Sf ( Fig.2b), this corresponds 
to C-CDW to IC-CDW to superfluid transition shown in 
Fig. la. Transferring back to the original BLQH problem, 
the small imbalance will first drive the C-PSDW to the 
In-commensurate pseudo-spin density wave (IC-PSDW), 
then drive a 1st order transition from the IC-PSDW to 
the ESF shown in Fig. lb. 



2Tr5f5, 



a=a/b 

where A^j, is a non-compact U{1) gauge field. Upto the 
quartic level, Eqn. 5 is the same in square lattice and in 
honeycomb lattice. If r > 0, the system is in the super- 
fluid state < ipi >= for every I = a/b. If r < 0, the 
system is in the insulating state < tjji >7^ for at least 
one I. We assume r < in Eqn. 5, so the system is in the 
insulating state. 71 > ( 71 < ) corresponds to the 
Ising ( or Easy-plane ) limit. The insulating state takes 
the CDW state ( or valence bond sohd (VBS) state ). 

In the balanced case Sf = 0, the SF to the VBS tran- 
sition in the easy plane limit was argued to be 2nd or- 
der through a novel deconfined quantum critical point 
[10]. However, the boson Hubbard model Eqn.4 on the 
PSDW side corresponds to the Ising limit in the dual 
model Eqn.5, therefore 71 > 0. The SF to the CDW 
transition in the Ising limit is first order. This is consis- 
tent with the first order ESF to PSDW transition driven 
by the collapsing of magnetoroton minimum studied in 
[9]. 

In the CDW order side, the mean field solution is 
ipa = l,fpb = 01 vice versa. In the slightly imbalance 
case Sf 7^ 0, setting = in Eqn.5 leads to: 

£ = 1(9^ - a^)Val' + rlVal' + 70lV'a|^ + " " " 

+ \ie^,.xd.Ax-2wSf5^rf (6) 
















ESF 


/iC-t'SDW 


FL 




./-^c5sDW 





r<0, Y,>0 



(a) 



(b) 



Fig.l: (a) The phase diagram of the boson Hubbard model 
Eqn.4 slightly away from the half filling, ij, is the chemical po- 
tential, (b) The bias voltage V versus distance d phase diagram 
at zero temperature. IC-PSDW stands for the incommensurate 
PSDW. The dashed line is the experimental path investigated in 
[13]. All the transitions are first order transitions. 

As shown in Fig. lb, the bias voltage increases, the im- 
balance will first introduce interstitials in the top layer 
and vacancies in the bottom layer, namely, turn the C- 
PSDW into the IC-PSDW whose charge distributions arc 
shown in Fig. 2a and 2b respectively, then the whole IC- 
PSDW melts into the ESF through a 1st order transition. 




Fig. 2: (a) The charge distribution of the PSDW in a square 
lattice. The dashed line stand for "up" pseudo-spins which take 
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sublatticc A, while "down" pseudo-spins take sublattice B. (b) The 
charge distribution of the IC-PSDW in a square lattice. The dot- 
ted line stand for "inverted" pseudo-spins which play the role of 
intcrstitials in the top layer or vacancies in the bottom layer. The 
number of interstitials in the top layer is equal to that of vacancies 
in the bottom layer. 

Disorders may smear all the 1st order transitions in 
Fig. lb into 2nd order transitions. The disorders may 
also transfer the long range lattice orders of the C-PSDW 
and IC-PSDW into short range ones. The fact make the 
observation of the commensurate and incommensurate 
lattice structures by light scattering experiment [14] dif- 
ficult. 

The dashed line in Fig. lb was investigated in a recent 
experiment [13]. But the first phase transition in Fig. lb 
was not paid attention in the experiment where the phase 
diagram was drawn against fixed charge imbalance in- 
stead of fixed bias voltage V. So the C-PSDW phase was 
crushed into the horizontal axis. A simple mean field ar- 
gument leads to the linear scaling of the second transition 
line y ~ d — del- A parabolic behavior /i^ ~ d — dd was 
found for the shape of the second transition at very small 
imbalances. We expect the disorders may transform the 
linear behavior to the parabolic one. In the presence 
of disorders, all the properties of the C-PSDW and IC- 
PSDW are consistent with the experimental observations 
in [13] on the intermediate phase at small imbalances. 

When the distance of the two layers is further increased 
to larger than a second critical distance dc2, then all the 
signature of the inter layer coherent state are lost, the 
two layers are decoupled into two separate v = 1/2 CF 
Fermi liquid state ( Fig. lb ). We expect that there is 
a level crossing and associated first order transition at 
dc2- When d > dc2, increasing the bias voltage may not 
transform the two decoupled FL state back into the ESF. 

(4) Discussions: By projecting into the LLL and then 
performing the HF approximation ( LLL-I-HF approach 
), the authors in [12] found that the transition at hz = 
0,d = del is an instability through a 1st order transi- 
tion to a pseudospin density wave state driven by the 
gap closing of magneto-roton minimum at a finite wave- 
vector. Starting from the ESF side, their numerical re- 
sults indicated that the imbalance increases the spin stiff- 
ness and also the critical distance dd ■ It is not known if 
the LLL-I-HF calculations are accurate in describing the 
transition from the ESF to the pseudo-spin density wave 
state. Our effective theory circninivcnts the difficulty as- 
sociated with the still not precisely known wavefunction 
at any finite d [15,16]. We study the effects of small 
imbalance from the PSDW side and map its effect as a 
chemical potential of hard core bosons with long range 



interactions hopping on a square lattice near half fill- 
ing, namely, mapping Fig. la to Fig. lb. Our upper phase 
boundary in Fig. lb is consistent with that achieved in 
[12] from the ESF side by the LLL+HF calculation. We 
also worked out the lattice structures of the C-PSDW 
and IC-PSDW and the whole physical picture along the 
dashed line in Fig. lb. It is not known how to apply the 
LLL-fHF theory in [12] to study the PSDW side. The 
efii'ective theory presented in this paper is complementary 
to and goes well beyond the previous LLL-I-HF calcula- 
tions. 
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Kong University in writing this paper. 
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